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Abstract 

Given a simplicial complex with weights on its simplices, and a nontrivial cycle on it, we are inter- 
ested in finding the cycle with minimal weight which is homologous to the given one. Assuming that the 
homology is defined with integer (Z) coefficients, we show the following (Theorem 5.2): 

For a finite simplicial complex K of dimension greater than p, the boundary matrix [dp+i] is 
totally unimodular if and only ifHp{L, Lq) is torsion-free, for all pure subcomplexes Lo, L 
in K of dimensions p and p+1 respectively, where Lq c L. 

Because of the total unimodularity of the boundary matrix, we can solve the optimization problem, which 
is inherently an integer programming problem, as a linear program and obtain an integer solution. Thus 
the problem of finding optimal cycles in a given homology class can be solved in polynomial time. This 
result is surprising in the backdrop of a recent result which says that the problem is NP-hard under Z2 
coefficients which, being ?i field, is in general easier to deal with. Our result implies, among other things, 
that one can compute in polynomial time an optimal [d — l)-cycle in a given homology class for any 
triangulation of an orientable compact d-manifold or for any finite simplicial complex embedded in M''. 
Our optimization approach can also be used for various related problems, such as finding an optimal 
chain homologous to a given one when these are not cycles. Our result can also be viewed as providing 
a topological characterization of total unimodularity. 



1 Introduction 

Topological cycles in shapes embody their important features. As a result they find apphcations in scientific 

studies and engineering developments. A version of the problem that often appears in practice is that given 
a cycle in the shape, compute the shortest cycle in the same topological class (homologous). For example, 
one may generate a set of cycles from a simplicial complex using the persistence algorithm [10] and then 
ask to tighten them while maintaining their homology classes. For two dimensional surfaces, this problem 
and its relatives have been widely studied in recent years; see, for example, [2, 3, 5, 6, 8]. A natural question 
is to consider higher dimensional spaces which allow higher dimensional cycles such as closed surfaces 
within a three dimensional topological space. High dimensional applications arise, for example, in the 
modeUng of sensor networks by Vietoris-Rips complexes of arbitrary dimension [7, 20]. Not surprisingly, 
these generalizations are hard to compute which is confirmed by a recent result of Chen and Freedman [4]. 
Notwithstanding this negative development, our result shows that optimal homologous cycles in any finite 
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dimension are polynomial time computable for a large class of shapes if homology is defined with integer 

coefficients. 

Let Khe a simplicial complex. Informally, a p-cycle in iiT is a collection of p-simplices whose bound- 
aries cancel mutually. One may assign a non-zero weight to each p-simplex in K which induces a weighted 
1-norm for each p-cycle in K. For example, the weight of a p-simplex could be its volume. Given any 
p-cycle c in K, our problem is to compute a p-cycle c* which has the minimal weighted 1-norm in the 
homology class of c. If some of the weights are zero the problem can still be posed and solved, except 
that one may not call it weighted 1-norm minimization. The homology classes are defined with respect to 
coefficients in an abelian group such as Q, M, Z, Z„ etc. Often, the group Z2 is used mainly because of 
simplicity and intuitive geometric interpretations. 

Chen and Freedman [4] show that under Z2 coefficients, computing an optimal p-cycle c* is NP-hard for 
p> 1. Moreover, their result implies that various relaxations may still be NP-hard. For example, computing 
a constant factor approximation of c* is NP-hard. Even if the rank of the p-dimensional homology group is 
constant, computing c* remains NP-hard for p > 2. The only settled positive case is a result of Chambers, 
Erickson, and Nayyeri [3] who show that computing optimal homologous loops for surfaces with constant 
genus is polynomial time solvable though they prove the problem is NP-hard if the genus is not a constant. 

The above negative results put a roadblock in computing optimal homologous cycles in high dimensions. 
Fortunately, our result shows that it is not so hopeless - if we switch to the coefficient group Z instead of Z2, 
the problem becomes polynomial time solvable for a fairly large class of spaces. This is a little surprising 
given that Z is not a field and so seems harder to deal with than Z2 in general. For example, Z2-valued 
chains form a vector space, but Z- valued chains do not. 

The problem of computing an optimal homologous cycle (or more generally, chain) can be cast as a lin- 
ear optimization problem. Consequently, the problem becomes polynomial time solvable if the homology 
group is defined over the reals, since it can be solved by linear programming. Indeed this is the approach 
taken by Tahbaz-Salehi and Jadbabaie [20]. However, in general the optimal cycle in that case may have 
fractional coefficients for its simplices, which may be awkward in certain applications. One advantage of 
using Z is that simplices appear with integral coefficients in the solution. On the other hand, the linear 
programming has to be replaced by integer programming in the case of Z. Thus, it is not immediately clear 
if the optimization problem is polynomial time solvable. One issue in accommodating Z as the coefficient 
ring is that integral coefficients other than 0, 1, or —1 do not have natural geometric meaning. Neverthe- 
less, our experiments suggest that optimal solutions in practice may contain coefficients only in {—1, 0, 1}. 
Furthermore, as we show later, we can put a constraint in our optimization to enforce the solution to have 
coefficients in {—1, 0, 1}. 

Our main observation is that the optimization problem that we formulate can be solved by linear pro- 
gramming under certain conditions, although it is inherently an integer programming problem. It is known 
that a linear program provides an integer solution if and only if the constraint matrix has a property called 
total unimodularity. A matrix is totally unimodular if and only if each of its square submatrices has a de- 
terminant of 0, 1, or —1. We give a precise topological characterization of the complexes for which the 
constraint matrix is totally unimodular. For this class of complexes the optimal cycle can be computed in 
time polynomial in the number of simpUces. Totally unimodular matrices have a well-known geometric 
characterization - that the corresponding constraint polyhedron is integral [15, Theorem 19.1]. Our result 
provides a topological characterization as well. 

We can allow several variations to our problem because of our optimization based approach. For ex- 
ample, we can probe into intermediate solutions; we can produce the chain that bounds the difference of 
the input and optimal cycles, and so forth. In fact, we can also find an optimal chain homologous to a 
given one when the chains are not cycles. In other words, we can leverage the flexibility of the optimiza- 
tion formulation by hnking results from two apparently different fields, optimization theory and algebraic 
topology. 
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2 Background 



Since our result bridges the two very different fields of algebraic topology and optimization, we recall some 
relevant basic concepts and definitions from these two fields. 

2.1 Basic definitions from algebraic topology 

Let K be a finite simplicial complex of dimension greater than p. A p-chain with Z coefficients in K is a 
formal sum of a set of oriented p-simplices in K where the sum is defined by addition in Z. Equivalently, it is 
an integer valued function on the oriented p-simplices, which changes sign when the orientation is reversed 
[14, page 37]. 

Two p-chains can be added by adding their values on corresponding p-simplices, resulting in a group 
Cp{K) called the p-chain group of K. The elementary chain basis for Cp{K) is the one consisting of 
integer valued functions that take the value 1 on a single oriented p-simplex, — 1 on the oppositely oriented 
simplex, and everywhere else. For an oriented p-simplex a, we use a to denote both the simplex and the 
corresponding elementary chain basis element. The group Cp{K) is free and abeUan. The boundary of an 
oriented p-simplex a = [vq, . . . ,Vp] is given by 

p 

9pO- = ^{-iy[vo, ..,Vi,..,Vp] , 

i=0 

where Vi denotes that the vertex Vi is to be deleted. This function on p-simplices extends uniquely [14, page 
28] to the boundary operator which is a homomorphism: 

dp-. Cp{K) Cp-i{K) . 

Like a linear operator between vector spaces, a homomorphism between free abelian groups has a unique 
matrix representation with respect to a choice of bases [14, page 55]. The matrix form of dp will be denoted 
[dp\. Let and {rjY-Z^ be the sets of oriented {p — 1)- and p-simplices respectively in K, ordered 

arbitrarily. Thus {crj} and {tj} also represent the elementary chain bases for Cp-\{K) and Cp{K) respec- 
tively. With respect to such bases [dp] is an m x n matrix with entries 0, 1 or — L The coefficients of dp tj 
in the Cp-i{K) basis become the column j (counting from 0) of [dp]. 

The kernel kei dp is called the group of p-cycles and denoted Zp{K). The image im5p+i forms the 
group of p-boundaries and denoted Bp{K). Both Zp{K) and Bp{K) are subgroups of Cp{K). Since 
dp o dp-^.l = 0, we have that Bp{K) C Zp{K), that is, all p-boundaries are p-cycles though the converse is 
not necessarily true. The p dimensional homology group is the quotient group Hp{K) = Zp{K) / Bp{K). 
Two p-chains c and c' in K are homologous if c = c' + Sp+i d for some (p + l)-chain d in K. In particular, 
if c = dp-^i d, we say c is homologous to zero. If a cycle c is not homologous to zero, we call it a non-trivial 
cycle. 

For a finite simplicial complex K, the groups of chains Cp{K), cycles Zp{K), and Hp{K) are all finitely 
generated abelian groups. By the fundamental theorem of finitely generated abelian groups [14, page 24] 
any such group G can be written as a direct sum of two groups G = F ®T where F = (Z © • • • © Z) and 
T = (Z/ti © • • • © Z/tfc) with ti > 1 and U dividing tj+i. The subgroup T is called the torsion of G. If 
T = 0, we say G is torsion-free. 

Let Lq be a subcomplex of a simplicial complex L. The quotient group Cp{L) / Cp{Lo) is called the 
group of relative chains of L modulo Lq and is denoted Cp{L, Lq). The boundary operator dp : Cp{L) — > 
Cp-i{L) and its restriction to Lq induce a homomorphism 

Cp{L,Lo)^Cp.i{L,Lo). 
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As before, we have dp^'^°^ od^^^^^ = 0. Writing Zp{L,Lq) = keTdp^'^°^ for relative cycles and 

Bp[L,Lq) = imdp^'i^^ for relative boundaries, we obtain the relative homology group Hp{L,Lq) = 
Zp{L, Lo)/Bp{L, Lq). Sometimes, to distinguish it from relative homology, the usual homology Hp{L) is 
called the absolute homology group of L. 

2.2 Total unimodularity and optimization 

Recall that a matrix is totally unimodular if the determinant of each square submatrix is 0, 1, or —1. The 
significance of total unimodularity in our setting is due to the following result: 

Theorem 2.1. [24] Let A be an m x n totally unimodular matrix and b an integral vector, i.e., b G Z"*. 

Then the polyhedron P := {x G R" | Ax. = b, x > 0} is integral, meaning that V is the convex hull of 
the integral vectors contained in V. In particular, the extreme points (vertices) ofV are integral. Similarly 
the polyhedron Q := {x G | j4x > b} is integral. 

The following corollary shows why the above result is significant for optimization problems. Consider 
an integral vector b G Z"* and a real vector of cost coefficients f G M". Consider the integer linear program 

min f^x subject to ^x = b, x > and x G Z" . (1) 

Corollary 2.2. Let Abe a totally unimodular matrix. Then the integer linear program (1) can be solved in 
time polynomial in the dimensions of A. 

Proof. Relax the integer linear program (1) to a linear program by removing the integrality constraint x G 
Z". Then an interior point method for solving linear programs will find a real solution x* in polynomial 
time [15] if it exists, and indicates the unboundedness or infeasibility of the linear program otherwise. In 
fact, since the matrix A has entries 0, 1 or —1, one can solve the linear program in strongly polynomial time 
[21, 22]. That is, the number of arithmetic operations do not depend on b and f and solely depends on the 
dimension of A. One still needs to show that the solution x* is integral. 

If the solution is unique then it lies at a vertex of the polyhedron V and thus it will be integral because 
of Theorem 2.1. If the optimal solution set is a face of V which is not a vertex then an interior point method 
may at first find a non-integral solution. However, by [1, Corollary 2.2] the polyhedron V must have at least 
one vertex. Then, by [1, Theorem 2.8] if the optimal cost is finite, there exists a vertex of V where that 
optimal cost is achieved. Following the procedure described in [12], starting from the possibly non-integral 
solution obtained by an interior point method one can find such an integral optimal solution at a vertex in 
polynomial time. □ 

3 Problem formulation 

Let be a finite simplicial complex of dimension p or more. Given an integer valued p-chain x = 
X^™Q^ x{ai) Gi we use x G Z"* to denote the vector formed by the coefficients x{ai). Thus, x is the 
representation of the chain x in the elementary p-chain basis, and we will use x and x interchangeably. For 
a vector v G W^, the 1-norm (or £^-norm) ||v||i is defined to be \vi\. Let W be any real mxm diagonal 
matrix with diagonal entries Wi . Then, the 1 -norm of v, that is, 1 1 v 1 1 1 is J2- \ wi\\vi\. (If 1^ is a general 
mxm nonsingular matrix then || v|| i is called the weighted 1-norm of v.) The norm or weighted norm 
of an integral vector v G Z"* is defined by considering v to be in W^. We now state in words the problem 
of optimal homologous chains and later formalize it in (2): 

Given a p-chain cin K and a diagonal matrix W of appropriate dimension, the optimal ho- 
mologous chain problem (OHCP) is to find a chain c* which has the minimal 1-norm ||H^c*||i 
among all chains homologous to c. 
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Remark 3.1. In the natural case where simplices are weighted and the optimality of the chains is to be 
determined with respect to these weights, we may take W to be diagonal with Wi being the weight of 
simplex cTj. In our formulation some of the weights can be 0. Notice that the signs of the simplex weights 
are ignored in our formulation since we only work with norms. 

Remark 3.2. In Section 1 we surveyed the computational topology literature on the problem of finding 
optimal homologous cycles. The flexibihty of our formulation allows us to solve the more general, opti- 
mal homologous chain problem, with the cycle case being a special case requiring no modification in the 
equations, algorithm, or theorems. 

Remark 3.3. The choice of 1-norm is important. At first, it might seem easier to pose OHCP using 2- 
norm. Then, calculus can be used to pose the minimization as a stationary point problem when OHCP is 

formulated with only equality constraints which appear in (2) below. This case can be solved as a linear 
system of equations. By using 1-norm instead of 2-norm, we have to solve a linear program (as we will 
show below) instead of a linear system. But in return, we are able to get integer valued solutions when the 
appropriate conditions are satisfied. 

The formulation of OHCP is the weighted -optimization of homologous chains. This is very general 
and allows for different types of optimality to be achieved by choosing different weight matrices. For 
example, assume that the simphcial complex K of dimension greater than p is embedded in R'^, where 
d > p + 1. Let 1^ be a diagonal matrix with the i-th diagonal entry being the Euclidean p-dimensional 
volume of a p-simplex. This specializes the problem to the Euclidean £^ -optimization problem. The resulting 
optimal chain has the smallest p-dimensional volume amongst all chains homologous to the given one. If W 
is taken to be the identity matrix, with appropriate additional conditions to the above formulation, one can 
solve the ^-optimization problem. The resulting optimal solution has the smallest number of p-simplices 
amongst all chains homologous to c, as we show in Section 3.2. 

The central idea of this paper consists of the following steps: (i) write OHCP as an integer program 
involving 1-norm minimization, subject to linear constraints; (ii) convert the integer program into an integer 
linear program by converting the 1-norm cost function to a linear one using the standard technique of 
introducing some extra variables and constraints; (iii) find the conditions under which the constraint matrix 
of the integer linear program is totally unimodular; and (iv) for this class of problems, relax the integer linear 
program to a linear program by dropping the constraint that the variables be integral. The resulting optimal 
chain obtained by solving the linear program will be an integer valued chain homologous to the given chain. 



3.1 Optimal homologous chains and linear programming 

Now we formally pose OHCP as an optimization problem. After showing existence of solutions we refor- 
mulate the optimization problem as an integer linear program and eventually as a linear program. 

Assume that the number of p- and {p + l)-simplices in is m and n respectively, and let 1^ be a 
diagonal m x m matrix. Given an integer valued p-chain c the optimal homologous chain problem is to 
solve: 

' ' (2) 



min||W^x||i such that x = c -|- [5p+i] y, and x e Z"*, y G Z" . 
x,y 



In the problem formulation (2) we have given no indication of the algorithm that will be used to solve the 
problem. Before we develop the computational side, it is important to show that a solution to this problem 
always exists. 

Claim 3.4. For any given p-chain c and any matrix W, the solution to problem (2) exists. 
Proof. Define the set 

Uo := {||l^x||i I X = c + [dp+i] y, X G and y G Z"} . 



5 



We show that this set has a minimum which is contained in the set. Consider the subset C Uc defined by 

U'^ = {\\W^\\i I ||VFx||i < ||W"c||i, x = c+[5p+i]y, xGZ^'andy gZ"}. 
This set is finite since x is integral. Therefore, inf Uc = inf = min U^. □ 

In the rest of this paper we assume that is a diagonal matrix obtained from weights on simplices 
as follows. Let it; be a real- valued weight function on the oriented p-simplices of K and let W be the 
corresponding diagonal matrix (the i-th diagonal entry of W is w{ai) = Wi). 

The resulting objective function ||VFx||i = Yli l""^?! I^^il in (2) is not linear in Xi because it uses the 
absolute value of Xj. It is however, piecewise-linear in these variables. As a result, (2) can be reformulated 
as an integer linear program in the following standard way [1, page 18]: 

min ^ \wi \ {xf + xj) 

i 

subject to x+ — x^ = c + y (3) 

x+, x" > 
x+, X- G Z", y G Z". 

Comparing the above formulation to the standard form integer linear program in (1), note that the vector 
X in (1) corresponds to [x+, x~, y]^ in (3) above. Thus the minimization is over x+, x~ and y, and the 
coefficients of x'l and x~ in the objective function are \wi\, but the coefficients corresponding to yj are zero. 
The linear programming relaxation of this formulation just removes the constraints about the variables being 
integral. The resulting hnear program is: 

min ^ l^il {xf + x~) 

i 

subject to x"*" — x~ = c + y (4) 

x+, X" > 0. 

To use the result about standard form polyhedron in Theorem 2.1 we can eliminate the free (unrestricted 
in sign) variables y by replacing these by y+ — y~ and imposing the non-negativity constraints on the 
new variables [I, page 5]. The resulting linear program has the same objective function, and the equality 
constraints: 

X+-X- =c+[5p+i] (y+-y-), (5) 

and thus the equahty constraint matrix is [l — / —B B\, where B = We now prove a result 

about the total unimodularity of this matrix. 

Lemma 3.5. If B = [dp+i] is totally unimodular then so is the matrix \l —I —B S] . 

Proof. The proof uses operations that preserve the total unimodularity of a matrix. These are fisted in [15, 
page 280]. If 5 is totally unimodular then so is the matrix \—B S] since scalar multiples of columns of 
B are being appended on the left to get this matrix. The full matrix in question can be obtained from this 
one by appending columns with a single ±1 on the left, which proves the result. □ 

As a result of Corollary 2.2 and Lemma 3.5, we have the following algorithmic result. 

Theorem 3.6. If the boundary matrix [9^+1 ] of a finite simplicial complex of dimension greater than p is 
totally unimodular, the optimal homologous chain problem (2) for p-chains can be solved in polynomial 
time. 
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Proof. We have seen above that a reformulation of OHCP (2), without the integrality constraints, leads 
to the linear program (4). By Lemma 3.5, the equality constraint matrix of this linear program is totally 
unimodular. Then by Corollary 2.2 the linear program (4) can be solved in polynomial time, while achieving 
an integral solution. □ 

Remark 3.7. One may wonder why Theorem 3.6 does not work when Z2-valued chains are considered 
instead of integer-valued chains. We could simulate Z2 arithmetic while using integers or reals by modify- 
ing (2) as follows: 

min||VFx||i such that x + 2 u = c + [Sp+i] y, and x G {0, 1}'", u e Z"*, y G . (6) 
x,y 

The trouble is that the coefficient 2 of u destroys the total unimodularity of the constraint matrix in the linear 
programming relaxation of the above formulation, even when [Sp+i] is totally unimodular. Thus we cannot 
solve the above integer program as a linear program and still get integer solutions. 

Remark 3.8. We can associate weights with {p + l)-simplices while formulating the optimization prob- 
lem (2). Then, we could minimize ||W^z||i where z = [x, y]-^. In that case, we obtain a p-chain c* ho- 
mologous to the given chain c and also a {p + 1) -chain d whose boundary is c* — c and the weights of c* 
and d together are the smallest. If the given cycle c is null homologous, the optimal y would be an optimal 
{p + 1) -chain bounded by c. 

Remark 3.9. The simplex method and its variants search only the basic feasible solutions (vertices of the 
constraint polyhedron), while choosing ones that never make the objective function worse. Thus if the 
polyhedron is integral, one could stop the simplex method at any step before reaching optimality and still 
obtain an integer valued homologous chain whose norm is no worse than that of the given chain. 



3.2 Minimizing the number of simplices 

The general weighted £^ -optimization problem (2) can be specialized by choosing different weight matrices. 
One can also solve variations of the OHCP problem by adding other constraints which do not destroy the 
total unimodularity of the constraint matrix. We consider one such specialization here - that of finding a 
homologous chain with the smallest number of simplices. 

If the matrix W is chosen to be the identity matrix, then one can solve the -optimization problem 
by solving a modified version of the ^^-optimization problem (2). One just imposes the extra condition 
that every entry of c and x be in {—1, 0, 1}. With this choice of Vl^ = / and with c G {—1, 0, 1}"*, the 
problem (2) becomes: 

min ||x||i such that x = c -|- y, andx G {-1,0, l}"", y G Z" . (7) 

Theorem 3.10. For any given p-chain c G { — 1,0, 1}™, a solution to problem (7) exists. Furthermore, 
amongst all x homologous to c, the optimal homologous chain x* has the smallest number of nonzero 
entries, that is, it is the £^ -optimal homologous chain. 

Proof. The proof of existence is identical to the proof of Claim 3.4. The condition that c takes values in — 1, 
0, 1 ensures that at least x = c can be taken as the solution if no other homologous chain exists. For the 

£° -optimality, note that since the entries of the optimal solution x* are constrained to be in {—1, 0, 1}, the 
1-norm measures the number of nonzero entries. Thus the 1-norm optimal solution is also the one with the 
smallest number of non-zero entries. □ 

Remark 3.11. Note that even with the given chain c taking values in {—1, 0, 1}, without the extra constraint 
that X G {—1, 0, 1}™ (rather than just x G Z™), the optimal 1-norm solution components may take values 
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outside { — 1, 0, 1}. For example, consider the simplicial complex K triangulating a cylinder which is shaped 
like an hourglass. Let ci and C2 be the two boundary cycles of the hour glass so that ci + C2 is not trivial. Let 
z be the smallest cycle around the middle of the hour glass which is homologous to each of ci and C2. Since 
c\ + C2 = 2z, the optimal cycle homologous to c\ + C2 has values 2 or —2 for some edges even if ci and C2 
have values only in {—1, 0, 1} for all edges. It may or may not be true that the number of nonzero entries 
is minimal in such an optimal solution. We have not proved it either way. But Theorem 3.10 provides a 
guarantee for computing ^^-optimal solution when the additional constraints are placed on x. 

The linear programming relaxation of problem (7) is 

min ^ {xj + xr) 

i 

subject to x''" — x~ = c + y (8) 

x+, x~ < 1 
x+, x~ > 0. 

One can show the integrality of the feasible set polyhedron by using slack variables to convert the inequal- 
ities x"*" < 1 and x~ < 1 to equalities and then using the V form of the polyhedron from Theorem 2.1. 
Equivalently, all the constraints can be written as inequalities and the Q polyhedron can be used. For a 
change we choose the latter method here. Writing the constraints as inequalities, in matrix form the con- 
straints are 
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(9) 



where B = [5p+i]. Then analogously to Lemma 3.5 and Theorem 3.6 the following are true. 

Lemma 3.12. IfB = is totally unimodular then so is the constraint matrix in (9). 

Theorem 3.13. If the boundary matrix [5p+i] of a finite simplicial complex of dimension greater than p 
is totally unimodular, then given a p-chain that takes values in { — 1, 0, 1}, a homologous p-chain with the 
smallest number of non-zeros taking values in {—1, 0, 1} can be found in polynomial time. 

In subsequent sections, we characterize the simplicial complexes for which the boundary matrix 
is totally unimodular. These are the main theoretical results of this paper, formahzed as Theorems 4.1, 5.2, 
and 5.7. 



4 Manifolds 

Our results in Section 5.1 are valid for any finite simplicial complex. But first we consider a simpler case - 
simplicial complexes that are triangulations of manifolds. We show that for finite triangulations of compact 
p-dimensional orientable manifolds, the top non-trivial boundary matrix [d^ is totally unimodular irrespec- 
tive of the orientations of its simplices. We also give examples of non-orientable manifolds where total 
unimodularity does not hold. Further examination of why total unimodularity does not hold in these cases 
leads to our main results in Theorems 5.2. 
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4.1 Orientable manifolds 



Let K he a. finite simplicial complex that triangulates a {p + 1) -dimensional compact orientable manifold 
M. As before, let [dp+i] be the matrix corresponding to dp+i : Cp+i{K) Cp{K) in the elementary chain 
bases. 

Theorem 4.1. For a finite simplicial complex triangulating a {p + 1)- dimensional compact orientable man- 
ifold, [dpj^i] is totally unimodular irrespective of the orientations of the simplices. 

Proof. First, we prove the theorem assuming that the [p + 1) -dimensional simplices of K are oriented 
consistently. Then, we argue that the result still holds when orientations are arbitrary. 

Consistent orientation of 1) -simplices means that they are oriented in such a way that for the {p+ 1)- 
chain c, which takes the value 1 on each oriented {p + l)-simplex in K, dp+i c is carried by the topological 
boundary d M of M. If M has no boundary then Sp+i c is 0. It is known that consistent orientation 
of [p + 1) -simplices always exists for a finite triangulation of a compact orientable manifold. Therefore, 
assume that the given triangulation has consistent orientation for the {p + l)-simpUces. The orientation of 
the p- and lower dimensional simplices can be chosen arbitrarily. 

Each p-face r is the face of either one or two (p + 1) -simplices (depending on whether r is a boundary 
face or not). Thus the row of [dp+i] corresponding to r contains one or two nonzeros. Such a nonzero entry 
is 1 if the orientation of r agrees with that of the corresponding {p + l)-simplex and —1 if it does not. 

Heller and Tompkins [13] gave a sufficient condition for the unimodularity of { — 1 , , 1 }-matrices whose 
columns have no more than two nonzero entries. Such a matrix is totally unimodular if its rows can be di- 
vided into two partitions (one possibly empty) with the following condition. If two nonzeros in a column 
belong to the same partition, they must be of opposite signs, otherwise they must be in different row par- 
titions. Consider the transpose of [9p+i]. Each column of [Sp+i]-^ contains at most two nonzero 
entries, and if there are two then they are of opposite signs because of the consistent orientations of the 
{p + 1) -dimensional simphces. In this case, the simple division of rows into two partitions with one con- 
taining all rows and the other empty works. Thus and hence [Op+i] is totally unimodular. 

Now, reversing the orientation of a (p -|- l)-simplex means that the corresponding column of [9p+i] 
be multiplied by —1. This column operation preserves the total unimodularity of Since any arbi- 

trary orientation of the {p + 1) -simplices can be obtained by preserving or reversing their orientations in a 
consistent orientation, we have the result as claimed. □ 

As a result of the above theorem and Theorem 3.6 we have the following result. 

Corollary 4.2. For a finite simplicial complex triangulating a (p + 1) -dimensional compact orientable 
manifold, the optimal homologous chain problem can be solved for p-dimensional chains in polynomial 
time. 

The result in Corollary 4.2 when specialized to W'^^ also appears in [19] though the reasoning is differ- 
ent. 

4.2 Non-orientable manifolds 

For non-orientable manifolds we give two examples which show that total unimodularity may not hold in 
this case. We also discuss the role of torsion in these examples in preparation for Theorem 5.2. 

Our first example is the Mobius strip and the second one is the projective plane. SimpUcial complexes 
for these two non-orientable surfaces are shown in Figure 1. The boundary matrices [^2] for these simplicial 

complexes are given in the Appendix in (1 1) and (12). 

Let M be the Mobius strip. We consider its absolute homology Hi{M) and its relative homology 
Hi{M, d M) relative to its boundary. Consult [14, page 135] to see how the various homology groups are 
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Figure 1: Triangulations of two non-orientable manifolds, shown as abstract simplicial complexes. The left figure 
shows a triangulation of the Mobius strip and the right one shows the projective plane. The numbers are the edge and 
triangles numbers. These correspond to the row and column numbers of the matrices (11) and (12). 



calculated using an exact sequence. We note that Hi (M) = Z, that is, its Hi group has no torsion. This can 
be seen by reducing the matrix (1 1) in the Appendix to Smith normal form (SNF). The SNF for the matrix 
consists of a 6 X 6 identity matrix on the top and a zero block below, which implies the absence of torsion. 

Let K be the simplicial complex triangulating M. Consider a submatrix S of the matrix [82] shown in 
Appendix as (11). This submatrix is formed by selecting the columns in the order 5, 4, 3, 2, 1, 0. From 
the matrix thus formed, select the rows 0, 3, 8, 9, 10, 2 in that order. This selection of rows and columns 
corresponds to all the triangles and the edges encountered as one goes from left to right in the Mobius 
triangulation shown in Figure 1 . The resulting submatrix is 

1" 






1 

1 -1 

The determinant of this matrix is —2 and this shows that the boundary matrix is not totally unimodular. 
The SNF for this matrix, it turns out, does reveal the torsion. This matrix S is the relative boundary matrix 
gi^'^o) where L = K and Lq are the edges in d M. The SNF has I's along the diagonal and finally a 2. This 
is an example where there is no torsion in the absolute homology but some torsion in the relative homology 
and the boundary matrix is not totally unimodular. We formulate this condition precisely in Theorem 5.2. 

The matrix [^2] given in Appendix as (12) for the projective plane triangulation is much larger. But it is 
easy to find a submatrix with determinant greater than 1 . This can be done by finding the Mobius strip in the 
triangulation of the projective plane. For example if one traverses from top to bottom in the triangulation of 
the projective plane in Figure 1 the triangles encountered correspond to columns 6, 9, 3, 8, 4 of (12) and the 



S = 



1 











-1 


1 











-1 


1 











-1 


1 











-1 
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edges correspond to rows 5, 1 1, 13, 12, 7. The corresponding submatrix is 



-1 
-1 








1 
1 









1 
1 










1 
1 



1' 






1 



and its determinant is —2. Tiius the boundary matrix (12) is not totally unimodular. Again, we observe 
that there is relative torsion in Hi{L, Lq) for the subcomplexes corresponding to the selection of S from 
[82]. Here L consists of the triangles specified above, which form a Mobius strip in the projective plane. 
The subcomplex Lq consists of the edges forming the boundary of this strip. This connection between 
submatrices and relative homology is examined in the next section. 

5 Simplicial complexes 

Now we consider the more general case of simplicial complexes. Our result in Theorem 5.2 characterizes 
the total unimodularity of boundary matrices for arbitrary simplicial complexes. Since we do not use any 
conditions about the geometric realization or embedding in M" for the complex, the result is also valid for 
abstract simplicial complexes. As a corollary of the characterization we show that the OHCP can be solved 
in polynomial time as long as the input complex satisfies a torsion-related condition. 

5.1 Total unimodularity and relative torsion 

Let K he a. finite simplicial complex of dimension greater than p. We will need to refer to its subcomplexes 
formed by the union of some of its simplices of a specific dimension. This is formalized in the definition 
below. 

Definition 5.1. A pure simplicial complex of dimension p is a simplicial complex formed by a collection of 
^?-simplices and their proper faces. Similarly, a pure subcomplex is a subcomplex that is a pure simpUcial 
complex. 

An example of a pure simphcial complex of dimension p is one that triangulates a p-dimensional man- 
ifold. Another example, relevant to our discussion, is a subcomplex formed by a collection of some p- 
simplices of a simplicial complex and their lower dimensional faces. 

Let L C K be a pure subcomplex of dimension p + 1 and Lq C L be a pure subcomplex of dimension 
p. If [dp+i\ is the matrix representing dp+i : Cp+i(i<') — > Cp{K), then the matrix representing the relative 
boundary operator 



is obtained by first including the columns of corresponding to {p + l)-simplices in L and then, from 

the submatrix so obtained, excluding the rows corresponding to the p-simplices in Lq and any zero rows. 
The zero rows correspond to p-simplices that are not faces of any of the {p + l)-simphces of L. 

As before, let be the matrix of 9p+i in the elementary chain bases for K. Then the following 

holds. 

Theorem 5.2. is totally unimodular if and only ifHp{L, Lq) is torsion-free, for all pure subcomplexes 

Lq,L of K of dimensions p andp+ 1 respectively, where Lq C L. 



d'X'r-Cp+i{L,Lo)^Cp{L,Lo), 
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Proof. (=^>) We show that if Hp{L, Lq) has torsion for some L, Lq then [dp+i] is not totally unimodular. 

to Smith normal form using 



Let 



aiL,Lo) 



be the corresponding relative boundary matrix. Bring 
the reduction algorithm [14] [pages 55-57]. This is a block matrix 



D 6 


where D = diag((ii , . . . , d;) is a diagonal matrix and the block row or column of zero matrices shown above 
may be empty, depending on the dimension of the matrix. Recall that di are integers and di> 1. Moreover, 
since Hp{L, Lq) has torsion, dk > I for some 1 < k < I. Thus the product di . . . dk is greater than 1. By a 
result of Smith [17] quoted in [15, page 50], this product is the greatest common divisor of the determinants 



of all X A; square submatrices of 



5, 



{L,Lq) 



. But this implies that some square submatrix of 



5, 



P+i 



, and 



hence of has determinant magnitude greater than 1. Thus is not totally unimodular. 



(<^) Assume that [9p+i] is not totally unimodular. We will show that then there exist subcomplexes Lq 
and L of dimensions p and (p + 1) respectively, with Lq C L, such that Hp{L, Lq) has torsion. Let S be 
a square submatrix of [9p+i] such that |det(S')| > 1. Let L correspond to the columns of [dp+i] that are 
included in S and let Sj, be the submatrix of [dp+i] formed by these columns. This submatrix Bl may 
contain zero rows. Those zero rows (if any) correspond to p-simphces that do not occur as a face of any of 
the {p + l)-simplices in L. In order to form S from Bl, these zero rows can first be safely discarded to form 
a submatrix B'j^. This is because det(5') / and so these zero rows cannot occur in S. 

The rows in B'^ correspond to ]5-simplices that occur as a face of some {p + 1) -simplex in L. Let Lq 
correspond to rows of B'j^ which are excluded to form S. Now S is the matrix representation of the relative 
boundary matrix dp^'^°\ Reduce S to Smith normal form. The normal form is a square diagonal matrix. 
Since the elementary row and column operations preserve determinant magnitude, the determinant of the 
resulting diagonal matrix has magnitude greater than 1. Thus at least one of the diagonal entries in the 
normal form is greater than 1. But then by [14, page 61] Hp{L, Lq) has torsion. □ 

Remark 5.3. The characterization appears to be no easier to check than the definition of total unimodularity 
since it involves checking every L, Lq pair. However, it is also no harder to check than total unimodularity. 
This leads to the following result of possible interest in computational topology and matroid theory. 

Corollary 5.4. For a simplicial complex K of dimension greater than p, there is a polynomial time algo- 
rithm for answering the following question: Is Hp{L,Lq) torsion-free for all subcomplexes Lq and L of 
dimensions p and + 1) such that Lq G L? 

Proof. Seymour's decomposition theorem for totally unimodular matrices [16], [15, Theorem 19.6] yields a 
polynomial time algorithm for deciding if a matrix is totally unimodular or not [15, Theorem 20.3]. That 
algorithm applied on the boundary matrix [9p+i] proves the above assertion. □ 

Remark 5.5. Note that the naive algorithm for the above problem is clearly exponential. For every pair L, Lq 
one can use a polynomial time algorithm to find the Smith normal form. But the number of L, Lq pairs is 
exponential in the number of p and {p + l)-simphces of K. 

Remark 5.6. The same polynomial time algorithm answers the question : Does Hp{L, Lq) have torsion for 
some pair L, Lq ? 
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5.2 A special case 

In Section 4 we have seen the special case of compact orientable manifolds. We saw that the top dimensional 
boundary matrix of a finite triangulation of such a manifold is totally unimodular. Now we show another 
special case for which the boundary matrix is totally unimodular and hence OHCP is polynomial time 
solvable. This case occurs when we ask for optimal d-chains in a simplicial complex K which is embedded 
in M'^+^. In particular, OHCP can be solved by linear programming for 2-chains in 3-complexes embedded 
in M^. This follows from the following result: 

Theorem 5.7. Let K be a finite simplicial complex embedded in W^'^^. Then, H^{L, Lq) is torsion-free for 
all pure subcomplexes Lq and L of dimensions d and d+1 respectively, such that Lq C L. 

Proof. We consider the {d + l)-dimensional relative cohomology group H'^+\L,Lo) (See [14] for exam- 
ple). It follows from the Universal Coefficient Theorem for cohomology [14, Theorem 53.1] that 

H'^+\L,Lo) = Rom{Hd+i{L,Lo), Z) Ext(i/^(L, Lq), Z) 

where Hom is the group of all homomorphisms from H^+i^L, Lq) to Z and Ext is the group of all of 
extensions between Hd,{L, Lq) and Z. These definitions can be found in [14, Chapter 5 and 7]. The main 
observation is that if H^{L,Lq) has torsion, Ext(i?(i(^, ^o), ^) has torsion and hence H'^^^{L,Lq) has 
torsion. 

On the other hand, by Alexander Spanier duality [18, page 296] 

H''+\L,Lq) = Ho{R''+' \ |Lo|,M'^+^ \ \L\) 

where \L\ denotes the underlying space of L. Since 0-dimensional homology groups cannot have torsion, 
ffd+i cannot have torsion. We reach a contradiction. □ 

Corollary 5.8. Given a d-chain c in a weighted finite simplicial complex embedded in W^"^^, an optimal 
chain homologous to c can be computed by a linear program. 

Proof. Follows from Theorem 5.7, Theorem 5.2, and Theorem 2.2. □ 

5.3 Total unimodularity and Mobius complexes 

As another special case, we provide a characterization of the total unimodularity of {p+ 1) -boundary matrix 
of simplicial complexes in terms of a forbidden complex called Mobius complex, for p < 1. In contrast to 
the previous characterization (in terms of relative homology of K), we directly employ certain results on 
totally unimodular matrices to derive this characterization in terms of submatrices called cycle matrices. We 
show in Theorem 5.13 that the {p + 1) -boundary matrix of a finite simplicial complex for p < 1 is totally 
unimodular if and only if the input complex does not have a (p + 1) -dimensional Mobius complex as a 
subcomplex. In particular, this observation along with Theorem 5.2 implies that a 2-complex does not have 
relative torsion if and only if it does not have a Mobius complex in it. We also demonstrate by example that 
this result does not generalize to higher values of p. 

Definition 5.9. A{p + 1) -dimensional cycle complex is a sequence ao, . . . , CTk-i of {p + l)-simplices such 
that ai and aj have a common face if and only ifj = {i + l) (mod A;) and that common face is a p-simplex. 
Such a cycle complex triangulates a{p + 1) -manifold. We call ita{p + 1) -dimensional cylinder complex if 
it is orientable and a.{p + 1) -dimensional Mobius complex if it is nonorientable. 
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Definition 5.10. For k > 2, a k x k matrix C is called a k-cycle matrix (k-CM) if Cij G {—1, 0, 1}, and C 
has the following form up to row and column permutations and scalings by —1: 



C = 



10 
1 1 
1 1 









0/3' 





1 
1 1 

oil 



(10) 



A k-CM with /3 = (—1)'^ is termed a cylinder cycle matrix (fc-CCM), while one with j3 = (—1)'^+^ is 
termed a Mobius cycle matrix (A;-MCM). We will refer to the form shown in (10) as the normal form cycle 
matrix. 

As an example, consider a triangulation K of a Mobius strip with k > 5 triangles shown in Figure 2. Let 
Kq be the complex for the boundary of the Mobius strip. In the figure, Kq consists of the horizontal edges. 



Then the relative boundary matrix 



of the Mobius strip K modulo its edge Kq is a A;-MCM. The 

orientations of triangle Tk-i and that of the terminal edge eo are opposite if k is even, but the orientations 
agree if k is odd, giving /3 = (—1)'^+^. Note that in Section 4.2, the submatrix S of the boundary matrix of 
the Mobius strip was such a relative boundary matrix and it is an example of a 6-MCM. Another example in 
that section was the 5-MCM obtained from the boundary matrix of the projective plane. 

Similarly, we observe a A;-CCM as the relative boundary-2 matrix of a cylinder triangulated with k 
triangles, modulo the cylinder's edges. Reversing the orientation of an edge or a triangle results in scaling 
the corresponding row or column, respectively, of the boundary matrix by —1. These examples motivate the 
names "Mobius" and "cylinder" matrices - a cycle matrix can be interpreted as the relative boundary matrix 
of a Mobius or cylinder complex. So, we have the following result. 

Lemma 5.11. Let K be a finite simplicial complex of dimension greater than p. The boundary matrix [5p+i] 
has no fc-MCM/or any k >2 if and only ifK does not have any (p + l)-dimensional Mobius complex as a 
subcomplex. 





Figure 2: Triangulation of a Mobius strip with k triangles. 



It is now easy to see that the absence of Mobius complexes is a necessary condition for total unimodu- 
larity. We show that this condition is also sufficient for 2- or lower dimensional complexes. We first need 
the simple result that an MCM is not totally unimodular. 

Lemma 5.12. Let C bea k-CM. Then det C = if it is a k-CCM, and |det C\ = 2 if it is a k-MCM. 

Proof. The matrix C can always be brought into the normal form with a series of row and column exchanges 

and scalings by —1. Note that these operations preserve the value of |detC|. Now assume that C has 
been brought into the normal form and call that matrix C . We expand along the first row of C to get 
det C" = 1 + {-if +^13, and the claim follows. □ 
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Theorem 5.13. For p < 1, [dp+i] is totally unimodular if and only if the simplicial complex K has no 
Mobius subcomplex of dimension p + 1. 

Proof. (=>) If there is a Mobius subcomplex of dimension p + 1 in K, then by Lemma 5.11 an MCM 
appears as a submatrix of That MCM is a certificate for not being totally unimodular since 

its determinant has magnitude 2 by Lemma 5.12. 

(<^) Let K have no Mobius subcomplexes of dimension p+1. Then by Lemma 5. 1 1, there are no MCMs as 
submatrices of [Op+i]. Truemper [23, Theorem 28.3] has characterized all minimally non-totally unimodular 
matrices, i.e., matrices that are not totally unimodular, but their every proper submatrix is totally unimodular. 
These matrices belong to two classes, which Truemper denotes as Wi and W-j. MCMs constitute the first 
class #i. A minimally non-totally unimodular matrix W isinWr if and only if W has a row and a column 
containing at least four nonzeros each [23, Cor. 28.5]. Since p < 1, no column of can have four or 

more nonzeros, and hence no matrix from the class can appear as a submatrix. Hence [5p+i] is totally 
unimodular if K has no {p + 1) -dimensional Mobius subcomplexes. □ 

The necessary condition in Theorem 5.13 extends beyond 2-complexes as Remark 5.14 indicates. How- 
ever, we cannot extend the sufficiency condition; Remark 5.15 presents a counterexample. 

Remark 5.14. Note that the absence of Mobius subcomplexes is a necessary condition for to be totally 

unimodular for all p. More precisely, if the simplicial complex K of dimension greater than p has a Mobius 
subcomplex of dimension p+1 then is not totally unimodular. By Lemma 5.1 1, an MCM appears as 

a submatrix of [3p+i] in this case. Its determinant has magnitude 2 by Lemma 5.12, trivially certifying that 

[5p+i] is not totally unimodular. 

Remark 5.15. The characterization in Theorem 5.13 does not hold for higher values of p. We present a 3- 
complex which does not have a 3-dimensional Mobius subcomplex, but whose [83] is not totally unimodular. 
Consider the simplicial complex consisting of the following seven tetrahedra formed from seven points 
numbered 0-6: (0,1,2,3), (0,1,2,4), (0,1,2,5), (0,1,2,6), (0,1,3,4), (0,2,3,5), (1, 2, 3, 6). It can be 
verified that the 19 x 7 boundary matrix [^3] of this simplicial complex has the 7 x 7 matrix 



W = 



-1 



1 





-1 







-1 





-1 








1 





-1 




1 









-1 





1 





as a submatrix where dct(VF) = —2, certifying that [ds] is not totally unimodular. In fact, W is the only 
submatrix of [^3] which is not totally unimodular, and it belongs to the class #7 of minimally non-totally 
unimodular matrices. 



6 Experimental Results 

We have implemented our linear programming method to solve the optimal homologous chain problem. In 
Figure 3 we show some results of prehminary experiments. 

The top row in Figure 3 shows the computation of optimal homologous 1 -chains on the simplicial com- 
plex representation of a torus. The longer chain in each torus figure is the initial chain and the tighter shorter 
chain is the optimal homologous chain computed by our algorithm. The bottom row shows the result of the 
computation of an optimal 2-chain on a simplicial complex of dimension 3. The complex is the tetrahedral 
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triangulation of a solid annulus - a solid ball from which a smaller ball has been removed. Two cut-away 
views are shown. The outer surface of the sphere is the initial chain and the inner surface is computed as the 
optimal 2-chain. 

In these experiments we used the linear program (4). The initial chains used had values in {—1, 0, 1} on 
the simplices. In the torus examples for instance, the initial chain was everywhere except along the initial 
curve shown. The curve was given an arbitrary orientation and the values of the chain on the edges forming 
the curve were +1 or —1 depending on the edge orientation. In these examples, the resulting optimal chains 
were oriented curves, with values of ±1 on the edges along the curve. This is by no means guaranteed 
theoretically, as seen in the hour glass example in Remark 3.11. The only guarantee is that of integrality. 
However if it is essential that the optimal chain has values only in {—1, 0, 1} then the linear program (8) or 
it's Euclidean variant can be used, imposing the additional constraint on the values of the optimal solution 
X as shown in linear program (8). 




Figure 3: Some experimental results. See text for details on what is being computed here. 



7 Discussion 

Several questions crop up from our problem formulation and results. Instead of 1-norm ||l^x||i, we can 
consider minimizing ^ • Wi Xj. In this case, the weights appear with signs and solutions may be unbounded. 
Nevertheless, our result in Theorem 3.6 remains valid. Of course, in this case we do not need to introduce 
xf and since the objective function uses Xi rather than | Xj | . We may introduce more generalization in the 
OHCP formulation by considering a general matrix W instead of requiring it to be diagonal and then asking 
for minimizing ||iyx||i. We do not know if the corresponding optimization problem can be solved by a 
linear program. Can this optimization problem be solved in polynomial time for some interesting classes of 
complexes? 

We showed that OHCP under Z coefficients can be solved by linear programs for a large class of topo- 
logical spaces that have no relative torsion. This leaves a question for the cases when there is relative 
torsion. Is the problem NP-hard under such constraint? Taking the cue from our results, one can also ask 
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the following question. Even though we know that the problem is NP-hard under Z2 coefficients, is it true 
that OHCP in this case is polynomial time solvable at least for simplicial complexes that have no relative 
torsions (considered under Z)? The answer is negative since OHCP for surfaces in is NP-hard under Z2 
coefficients [3] even though they are known to be torsion-free. 

Even if the input complex has relative torsion, the constraint polyhedron of the linear program may still 
have vertices with integer coordinates. In that case, the linear program may still give an integer solution 
for chains that steer the optimization path toward such a vertex. In fact, we have observed experimentally 
that, for some 2-complexes with relative torsion, the linear program finds the integer solution for some input 
chains. It would be nice to characterize the class of chains for which the linear program still provides a 
solution even if the input complex has relative torsion. 

A related question that has also been investigated recently is the problem of computing an optimal 
homology basis from a given complex. Again, positive results have been found for low dimensional cases 
such as surfaces [11] and one dimensional homology for simplicial complexes [5, 9]. The result of Chen 
and Freedman [4] implies that even this problem is NP-hard for high dimensional cycles under Z2. What 
about Z? As in OHCP, would we have any luck here? 
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Appendix 



Boundary matrices for non-orientable surfaces 

The boundary matrices [82] for the Mobius strip and projective plane triangulations shown in Figure 1 are 
given below. The row numbers are edge numbers and the column numbers are triangle numbers which are 
displayed in Figure 1 . 











[^2] for Mobius strip 


















: 


1 : 


2 : 


3 : 


4 : 


5 















1 














1 










1 
















-1 













2 




-1 


1 






















3 
















1 


-1 










4 







-1 






















5 




1 

























6 



















1 










7 










-1 



















8 













1 


-1 













9 










1 


-1 
















10 







1 


-1 



















11 













1 




















[^2] for projective plane : 










: 


1 




2 : 


3 : 


4 : 


5 : 


6 : 


7 : 


8 : 


9 : 





-1 














-1 














1 





1 


1 























2 


1 




1 


























3 








-1 














1 








4 

















1 





-1 








5 














-1 





-1 











6 


-1 























1 





7. 














1 











-1 





8 

















-1 


1 











9 





1 























-1 


10 


: 





1 





-1 

















11 


: 

















-1 








1 


12 


: 








-1 














1 





13 


: 








1 

















-1 


14 


: 








-1 











1 









(11) 



(12) 
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